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Abundance analyses of helium-rich subluminous B stars 



o 

(N 
< 

co 

Oh 

6 

(N 
> 

in 
q 

o 
o 



X 



NaslimN. 1 *, C.S. Jeffery 1 !, A. Ahmad 1 , N.T. Behara 12 & T. §ahin 1 ' 3 

1 Armagh Observatory, College Hill, Armagh BT61 9DG, N. Ireland 

2 Institut dAstronomie et d'Astrophysique, Universite Libre de Bruxelles, Belgium 

3 'Department of Astronomy and The W.J. McDonald Observatory, The University of Texas at Austin, Austin, Texas 78712, USA 



Accepted Received ; in original form . 



ABSTRACT 

The connection between helium-rich hot subdwarfs of spectral types O and B (He-sdB) has 
been relatively unexplored since the latter were found in significant numbers in the 1980's. In 
order to explore this connection further, we have analysed the surface composition of six He- 
sdB stars, including LB 1766, LB 3229, SB 21 (= Ton-S 137 = BPS 29503-0009), BPS 22940- 
0009, BPS 29496-0010, and BPS 22956-0094. Opacity-sampled line-blanketed model at- 
mospheres have been used to derive atmospheric properties and elemental abundances. All 
the stars are moderately metal-poor compared with the Sun ([Fe/H]« —0.5). Four stars are 
nitrogen-rich, two of these are carbon-rich, and at least four appear to be neon-rich. The data 
are insufficient to rule out binarity in any of the sample. The surface composition and locus of 
the N-rich He-sdBs are currently best explained by the merger of two helium white dwarfs, or 
possibly by the merger of a helium white dwarf with a post-sdB white dwarf. C-rich He-sdBs 
require further investigation. 

Key words: stars: early-type, stars: subdwarfs, stars: chemically peculiar, stars: abundances, 
stars: evolution. 



1 INTRODUCTION 

Subdwarf B stars are low-mass core-helium burning stars with ex- 
tremely thin hydrogen envelopes. They behave as helium main- 
sequence stars of roughly half a solar mass. Their atmospheres 
are generally helium deficient; radiative levitation and gravitational 
settli ng combine to make helium sink below the hydrogen-rich sur- 
face dHeberll 1986b . 

However, almost 5% of the total subdw arf population com- 
prise stars with heliu m-rich atmospheres dGreen et alj 1 19861 : 
Ahmad & Jefferv 2006). The optical spectra of these stars are char- 
acterised by strong neutral helium lines and weak Hell lines; they 
exhibit a wide range of helium abundance and effective temper- 
atures (T c ff) similar to both sdB and sdO star s. They have been 
variously classified as sdOB, sdOC and sdOD dGreen et al. 1986) 
stars, but more recently as He-s dB and He-sdO stars dMoehler et al.l 
1990!; lAhmad& Jeffovll2004l) . The spectroscopic division con- 
cerns the relative strengths of Hel4471, Hell 4541 and H7 (itself 
a blend of H and Hell). Ro ughly speakin g, the division occurs for 
stars with T cB « 38 000 K lDrilling et"al] I2003I) . 

In general, He-sdB stars hav e lower surfac e gravities (g) than 
normal hydrogen- rich sdB stars dHeberl 120091) and have spectral 
characteristic s intermediate bet ween extreme heli um (EHe) stars 
dleffervll 19961) and He-sdO stars dNapiwotzkill2008l) . 

Most He-sdB stars show strong nitrogen (Nil and Nlll) lines 
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in their optical spe ctrum; these are referred to as N-rich by 
iDrilling et al.l d2003l) . A few also showing strong carbon (ClI and 
Cm) lines are labelled C-rich. 

The question posed by these stars is that of their evolutionary 
status. Do He-sdB and He-sdO stars form a single sequence? Why 
are there C-rich and N-rich stars? Why is there such a large range 
in hydrogen abundance? How are they related to other classes of 
evolved star, including normal s dB stars? Is the re a connection with 
any of the extreme helium stars dJeffervll2008a llbl)? 

Possible origins include: a late core flash of a single post- 
giant-branch helium star evolving toward the w hite dwarf sequence 
ILanz et alJl2004l: Mil ler Bertolami et alj|2008l): the merg er of two 
helium white dwarfs dlberj 1 1990l : ISaio & Jeffervl 12000): and the 
merg er of a helium white dwarf with a post-sdB star dlustham et al.l 
2010). All of these scenarios are likely to produce hot subdwarfs 
with He-rich and N-rich surfaces. The second has also been argued 
to lead to helium-poor "normal" sdB stars. All predict evolution 
tracks that commence with shell helium-ignition in a white dwarf. 
They take the star to a yellow-giant on a thermal time-scale and 
then to the helium main-sequence on a nuclear timescale. The de- 
tails of the tracks differ in respect of their initial conditions and the 
micro-physics adopted. Whether carbon is exposed is not clear. 

The goal is, if possible, to distinguish clearly the various types 
of observed He-sdB (and He-sdO) stars and to connect each to one 
of these diverse origins. 

The surface abundances of elements other than hydrogen and 
helium are therefore important indicators of previous evolution. But 
first it is necessary to establish in what ranges these abundances lie, 
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and hence to identify whether distinct groups exist. The challenge 
is that the numbers of He-sdB and He-sdO stars bright enough for 
fine analysis has hitherto been small, and very few, so far, have been 
found to be spectroscopically similar. The first of our studies was 
of PG 1544+488, which surprisingly turned out t o be a short-period 
binar y containing two helium-rich sdB stars jAhmad & J efferv 
|2004 . The second was for JL87 jAhmad et al.1 |2007|) . a rela- 
tively bright and only moderately helium-rich (nnc/nu ~ 0.4) 
and carbon-rich subdwarf . Far-ultraviole t spectra of t hese, and of 
LB 1766, were previously analyze d by Lanz et alj 1 2004). with 
quite different results to our own. lAhmad et alj 1 2007 ) demon- 
strated the importance of establishing T e s, log 9 and carbon abun- 
dances from optical spectra of Hel lines, i.e. relatively unblended 
lines with well-understood broadening theory, before extracting 
abundances of subordinate species. 

We are therefore systematically acquiring high-resolution high 
signal-to-noise optical spectroscopy of He-sdB stars. § 2 describes 
the observations used in this paper. These data are used to carry out 
detailed abundance analyses, making use of the latest generation of 
fully line-blanketed model atmospheres for appropriate mixtures 
(§,3). §4 presents the results for our programme stars, which are 
discussed in terms of the evolution models and analyses of related 
objects in § 5. 



2 OBSERVATIONS 

Spectra of several hydrogen-deficient stars were obtained with the 
University College London Echelle Spectrograph on the Anglo- 
Australian Telescope (A AT) in August 2005. These included a 
number of He-sdB stars, namely LB 1766, LB 3229, SB 21 (= Ton- 
S137 = BPS 29503-0009), BPS 22940-0009, BPS 29496-0010, 
and BPS 22956-0094. An extract from the observing log is shown 
in Table[T] which also indicate the signal-to-noise ratio of the com- 
bined spectra used here. UCLES was configured with the 31.6 lines 
mm -1 grating, the EEV2 detector, and a slit width of 1.09 mm. A 
central wavelength of 4340.02 A gives complete spectral coverage 
between 3820 and 5200 A, and a nominal resolution with this slit 
width R « 32 000. Exposures were broken into 1800 s segments in 
order to minimise cosmic-ray contamination. 

A preliminary analysis of three of these stars (LB 1766, 
SB 2 1, and BPS 22940 -0009) using the same data was given 
by jNaslimetal] |2QljJ). SB 2 1 w as originally identified a s 
an extremely helium -rich subdwarf ( Hunger & Kud ritzkj Il980l) : 
Hung eret alj J 198 lh found this star to be comparable with the 
helium-rich hot subdwarfs, CPD-31°1701 and TON-S 103. Abun- 
dances for LB 17 66 were previous ly obtained from a far-ultraviolet 
(FUV) spectrum jLanz et all 1 20041) ; experience has demonstrated 
that an ana lysis of the FUV sp ectrum alone can lead to system- 
atic errors jAhmad et al.1 120071) . Our AAT spectra show LB 1766 
and S B 21 to be nearly identical. Note that t he original analy- 
ses jLanz et alj2004lHunger & Kudritzkill98d) suggest these two 
stars to be quite different; thus a detailed contemporary comparison 
is important. 

BPS 22940-0009 and BPS 29496-0010 were identified as 
He-sdB stars and BPS 22956-0094 as an sdB star in the sur- 
Our spectra show BPS 22940-0009 



vey of iBeers et al.1 ([1992). 
and BPS 22956-0094 to be carbon-rich He-sdB stars with strong 
ClI, Cm, Nil, NlII, and Hel lines. In contrast, LB 1766, SB 21, 
BPS 29496-0010 and LB 3229 are carbon poor He-sdB stars, but 
with strong Nil, NlII, and Hel lines. 

The observations were reduced using a combination of 



Table 1. AAT/UCLES Observing Log 







a (2000) 5 


my 


u 1 (Mart) 


Seeing 


"CJijJ V / 


S/N 




BPS CS 22940-0009 




20 30 20 


-59 50 50 


13.8 


2005 08 26 09:25:40 


1.4" 


1800 


10 




2005 08 26 09:56:33 


1.4" 


1800 


10 




2005 08 26 10:27:29 


1.4" 


1800 


12 




2005 08 26 10:58:22 


1.4" 


1800 


10 




Mean 






21 




LB 1766 




04 59 19 


-53 52 52 


12.3 


2005 08 26 17:22:30 


2.0" 


1800 


35 




2005 08 26 17:53:23 


2.0" 


1800 


35 




2005 08 29 16:00:55 


1.7" 


1800 


20 




2005 08 29 16:31:49 


1.7" 


1800 


20 




Mean 






61 




BPS L.S 22956—0094 




22 16 56 


-64 31 51 




2005 Uo 2/ 1 1 AL. 1 y 


1.5" 


1800 


17 




or^c f\o t~i 11.11.11 

2005 Uo 27 12:13:13 


1.5" 


1800 


13 




2005 08 27 12:44:07 


1.5" 


1800 


20 




Mean 






28 




BPS CS 29496— 00 10 




23 34 02 


-28 51 38 


14.7 


2U0j 05 2/ 14.0/.JZ 


1.8" 


1800 


10 




2005 08 27 14:38:46 


1.8" 


1800 


10 




2005 08 27 15:09:40 


1.8" 


1800 


to 




2005 08 27 15:40:37 


1.5" 


1800 


to 




2005 08 27 16:21:17 


1.5" 


1800 


to 




2005 08 27 16:52:10 


1.5" 


1800 


10 




2005 08 27 17:23:04 


1.5" 


1800 


10 




2005 08 27 17:54:36 


1.5" 


1800 


to 




Mean 






25 




LB 3229 = JL 261 




01 47 17 


-51 33 39 


13.6 


2005 08 27 18:29:28 


1.5" 


1800 


17 




2005 08 27 19:00:23 


1.5" 


1800 


18 




Mean 






25 




SB 21 = PHL645 = Ton-S 137 


= BPS CS 29503-0009 








00 04 31 


-24 26 18 


13.9 


2005 08 28 14:46:45 


1.2" 


1800 


20 




2005 08 28 15:17:39 


1.2" 


1800 


20 




2005 08 28 15:48:33 


1.2" 


1800 


20 




2005 08 28 16:19:27 


1.2" 


1800 


19 




2005 08 28 16:50:21 


1.2" 


1800 


18 




2005 08 28 17:21:15 


1.2" 


1800 


18 




2005 08 28 17:52:09 


1.2" 


1800 


15 




2005 08 28 18:24:53 


1.4" 


1800 


17 




Mean 






52 





ECHOMOP routi nes jMills et alJuTOOrj) and bespoke echelle re- 
duction software iSahinll2008h ~ The sky-subtracted wavelength- 
calibrated spectrum was extracted to a 2-d format with each or- 
der represented by a single row. Continuum normalization was 
achieved by smoothing the 2-d spectrum to form a 2-d envelope 
spectrum, and then dividing by said envelope to remove most of 
the echelle blaze function. The smoothing procedure was adjusted 
to ensure that strong lines (Balmer or Hel 4471 for example) were 
avoided when producing the normalisation function. Order merging 
was carried out using the same 2-d envelope to provide the weights 
at each wavelength in each overlap interval. A final normalisation 
step was carried out in which the merged 1-d spectrum was divided 
by a low-order polynomial fitted to a set of continuum points de- 
fined manually. 

We found no radial velocity shifts amongst the repeat spectra 
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Figure 1. Theoretical model tits to Hel4471 in LB 3229 for model atmo- 
spheres with n He = 0.989, T cff = 40 000 K and log g = 4.0(0.5)6.0. 
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Figure 2. The loci of ionization equilibria for Nn/NlII and Hell, the profile 
fits to Hel lines, and the adopted solution for LB 3229. 



of individual stars, but only the data for LB 1766 were spread over 
a significant time interval (3 days). 

In order to obtain a spectrum with the highest possible signal- 
to-noise ratio, all of the individual spectra for each object were 
merged together, weighted appropriately for the signal-to-noise ra- 
tio in each spectrum. The combined spectra were velocity-shifted 
to the rest frame. 

In an iterative process, we successfully identified all of the sig- 
nificant absorption lines visible in the combined spectra. Where ob- 
servations allow (e.g. LB 1766 and SB 21), all known per mitted and 
forbidden lines of Hel can be id e ntified (cf. HP 14494 1: Underhill 
1 19661 : lHarrison & Jeffervl 1 19971 : iBeauchamp & Wesemaell 1 19981) . 
Hell4686 is present in all targets. Hydrogen Balmer lines are ev- 
ident in BPS 22956-0094. Hell4541 and other Hell Pickering lines 
are present in LB 3229 and BPS 29496-0010. In the remainder, 
an H/3/HeIl4859 blend is present; it is not obvious from H7 or 
Hell4541 which is dominant. In addition to Nil, III, the target spec- 
tra variously show lines due to ClI.III, Oil, Nell, Mgll, Aim, 
Sinn, IV and Sill. Identification charts are given in Figs. lAll - IA6l 
(online only). 



3 PHYSICAL PARAMETERS OF HE-SDBS 

The goal was to measure atmospheric physical parameters T e s, 
logg and elemental abundances for each star. We adopted two 
methods to determine T a g and log g. In the first we used ionisa- 
tion equilibria of prominent ions to determine T c h and the profiles 
of Stark-broadened Hel lines to determine log g. In the second we 
used the \ 2 -minimization package SFIT ( Ahmad & Jeff erv [2003b 
to determine T c a and log g simultaneously. 

3.1 Model atmospheres, line formation and spectral synthesis 

Grids of model atmospheres for hydrogen-deficient stars 
were calculated using the LTE line-blanketed code STERNE 
jBehara& Jeffervl l2006h which uses Opacity Project data for the 
continuous opacities, and treats line-blanketing through opacity 
sampling in a database of some 10 s atomic transitions. In this 
case, we adopted grids with 1/10 solar metallicity, relative helium 
abundances n Hc = 0.50, 0.70, 0.90, 0.95, 0.989 and 0.999 by 
numbeiQ, and assumed micro-turbulent velocities v t = 5 and 
lOkms" 1 . 

The choice of the 1/10 solar metallicity grid was adopted pri- 
marily because we were unable to obtain satisfactory fits for solar 
abundance models. The average abundances of Si and Mg, which 
are unlikely to have been affected by evolution, are sub-solar by 
« 0.5 dex (see below). 

Given a model atmosphere , the LTE radiative-transfer code 
SPECTRUM ( I Jeff erv et all200ll) may be used to compute a) indi- 
vidual line profiles and equivalent widths for given abundances and 
Vt, b) synthetic spectra for given wavelength ranges given the same 
information, and c) the abundances of ions from individual lines of 
given equivalent width and v t . 

For all of the programme stars except BPS CS 22956-0094, 
the Balmer lines were either weak or undetectable. Consequently 
we initially assumed a helium abundance nHc = 0.999 (nn = 0.0) 
and model atmospheres with v t = 5kms _1 as a starting approxi- 
mation. 

3.2 Micro-turbulent velocity 

Using approximate values for T c s and logg, we measured the 
micro-turbulent velocity from the equivalent widths of seven- 
teen Nil lines in LB 1766 and SB 21. Nitrogen abundances were 
calculated for micro-turbulent velocities in the range v t = 

(5) 20kms _1 using SPECTRUM. The micro-turbulent velocity 
determined by minimizing the scatter in the nitrogen abundance 
was v-t = 10 ± 4kms~ 1 . We used the measured value of v t in 
subsequent formal solution calculations for the ionisation equilib- 
rium and abundance measurements. For the instrumental profile, 
we adopted a Gaussian with FWHM = 2 resolution elements, cor- 
responding to 0.1 A or R « 45 000. Additional broadening was 
attributed to rotation broadening and measured as part of the SFIT 
solution. 

3.3 Ionisation equilibrium and Hel fitting 

Using model atmospheres with v t = lOkms" 1 , the ionisation 
equilibrium was established by balancing the abundances deter- 
mined from Nil and NlII lines in each star, as well as by fitting the 

1 equivalent to n H = 0.50, 0.30, 0.10, 0.05, 0.01 and 0.0, respectively. 
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equivalent width of Hell 4686. Values of T e g were determined for 
several fixed values of log g. Ideally, all ionisation equilibria and 
profile fits should converge at a single point in parameter space. In 
practice they rarely do, possibly because of departures from LTE, 
which become increasingly important at T e ff > 30 000 K espe- 
cially in certain strong Hel and Hell lines, possibly because of er- 
rors in equivalent width measurement, atomic data, or some other 
reason. 

Similarly, for a series of fixed values of T c g, we determined 
the value of log g by finding the best-fit theoretical profile for the 
Hel lines 4471, 4388 and 4922 A (Fig. [TJ. The non-diffuse line 
Hel 4121 is blended with nearby lines, and so was excluded from 
our analysis. 

The coincidence of Hel profile fits and the ionisation equilib- 
ria was used to determine the overall solution illustrated in Fig. [2] 
Since the Hell 4686 and Nil/Ill temperatures do not coincide, and 
since the systematics of these two diagnostics are not yet clear, we 
have chosen an unweighted mean to determine the ionization tem- 
perature. The error in this mean is given by the quadratic mean of 
the formal error in the Hell 4686 temperature and the standard de- 
viation in the Nil/Ill temperature. 

In practice, this procedure requires making a choice for the 
helium abundance. An estimate for who was obtained as described 
in the following section, and the grid helium abundance closest to 
that estimate was adopted. 

3.4 x 2 -minimization 

We measured the physical parameters T e g, logg and hydrogen 
abundance nn using the package SFIT which finds the best-fit solu- 
tion within a grid of synthetic spectra. Although 7iHc is a parameter 
of the model grid, at high helium abundances it is actually a proxy 
for riH = 1 — Ei>irii; the Hel lines are insensitive to abundance 
whilst, unusually for B stars, the Balmer lines vary strongly with 
abundance. SFIT works by fitting as many parameters and regions 
of spectrum simultaneously as the user chooses. Thus it measures 
T c s from the relative strengths of helium lines, the ionization equi- 
libria of all elements in the spectrum (e.g. Hel/Il, Nil/Ill), riHe (or 
71h) from the strengths of hydrogen and helium lines and log g from 
profiles of Stark-broadened lines. In our analysis, spectral regions 
blueward of 4050 A where the broad wings of Hel lines merge with 
one another, have been excluded because normalisation is difficult. 

The model grid used for the analysis of each star was 
a subset of a larger model grid in which T e g = 32 000 — 
40 000, 50, 000 K, log 3 = 4.00 - 6.00, v t = lOkms" 1 , and 
7i Hc = 0.949, 0.989, 0.999, except in the case of BPS 22956- 
0094, where n Hc = 0.500, 0.699, 0.899 was used. The grid spac- 
ings were 5T e f{ = 2000 K and Slogg = 0.5. Since no grid was 
available for 7iHc = 0.949 with v t = lOkms -1 , we replaced this 
with a grid having v t — 5 kms" . The spectral fitting was done it- 
eratively. In the initial iteration SFIT was run with three free param- 
eters; T c ff, logg and nHc- The value of nHc was noted and fixed. 
In the final iteration T e g and log g were solved simultaneously. 

3.5 Atmospheric parameters 

The atmospheric parameters of each star measured separately using 
SFIT and ionisation equilibrium are given in Table [2] Both sets of 
results agree with one another to within the formal errors. In the 
case of LB 1766, T e g is lower than in both previous studies, but 
close to that for SB 21, as would be indicated by the spectral sim- 
ilarity of these two stars. The model atmospheres are substantially 



improved, incorporating more appropriate line blanketing than in 
any previous studies, and this does acco unt for significant shifts in 
T e ff in hydrogen-deficient atmospheres (Behara & Jeffery 2006). 

3.6 Chemical abundances 

Having measured T e s, \ogg and nuc for each star using two dif- 
ferent methods, we chose the grid model atmosphere closest to 
these measured values (labelled "Model" in Table|2](. We measured 
equivalent widths of all C, N, O, Ne, Mg, Al, Si, and S lines for 
which we had atomic data. SPECTRUM can compute a curve of 
growth for any given spectral line; given an equivalent width it 
will then return the elemental abundances for that line. Table [ATI 
gives the adopted oscillator strengths (<?/), measured equivalent 
widths, and line abundances. Abundances are given in the form 
€i = log rii + c where log E^n^ = log E^rii© = 12.15 and a% 
are atomic weights. This form conserves values of t% for elements 
whose abundances do not change, even when the mean atomic mass 
of the mixture changes substantially. 

Mean abundances for each element are reported in Table[3j in 
general, the errors represent the standard deviation of the line abun- 
dances about the mean. However, for Mg, Al, and S, the error also 
includes the error on the equivalent width measurement estimated 
from the continuum noise. The errors in abundance Set due to a 
representative systematic change in <5T e g or 5 log g are shown, for 
three stars, in Table lA2l 

The hydrogen abundances adopted previously were further re- 
fined by starting with a model spectrum for each star defined by 
the best model (Table [2} and abundances (Table [3}, and by using 
SFIT to solve for the hydrogen abundance by fitting the Balmer 
lines only. In practice, only an upper limit of nn < 0.001 could 
be established for three stars, while the remaining three have the 
abundances shown in Table[3] 

Table[3]compares the elemental abundances thus derived with 
those previously published for LB 1766, with those of more H- 
rich "He-sdB" stars JL 87 and LS IV- 14° 116, two extreme helium 
stars V652 Her and BX Cir, and the Sun. These will be discussed in 
the next section. 

Theoretical spectra computed using the adopted grid atmo- 
sphere "Model" (TableO and the adopted mean abundances (Table 
O are shown over-plotted on the observed spectra in Figs. lAl I IA6I 

We note the following: 

1) Using silicon (five to eight lines), magnesium (one line), alu- 
minium (two lines) and sulphur (two - three lines) as proxies for 
overall metallicity, the group is metal poor by « 0.5 ± 0.2 dex 
compared with the Sun. We have not yet identified any iron lines in 
the optical spectra or analysed the FUSE spectrum of LB 1766. 

2) The majority are hydrogen-deficient. The Balmer lines may be 
blended with weak Hell lines. The hydrogen abundances are mea- 
sured by % 2 minimisation in the model grid; errors are estimated. 
The exception is BPS CS 22956-0094. 

3) All stars are nitrogen-rich (+0.46 ± 0.11 dex) compared with 
the Sun, and significantly so (~ +1 dex) after allowing for their 
low metallicity. 

4) After correcting for metallicity, the group ranges from very 
carbon-rich (BPS 22940-0009 and BPS 22956-0094: « +0.8 dex) 
to very carbon-poor (SB 21 and BPS CS 22496-0010: ss -1.2 
dex). 

5) In the two C-rich stars, the broad absoprtion at 4618A is not sat- 
isfactorily reproduced in the mode l. This an omalous line was first 
identified in H-d eficient spectra bylKlemolal dl96ll) . and discussed 
most recently bv lLeuenhagen et afl i 1996h . 
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Table 2. Atmospheric parameters 



Star 




T cff (K) 


log g 




(kms^ 1 ) 


v sin i 
(km s ) 


Source 


BPS CS 22940-0009 


33 700 ± 800 


4.7 ± 0.2 


0.993 




4 ± 3 


SFIT 






34 150 ± 1 700 


4.8 ± 0.2 








Ionization equilibrium 






34 000 


4.5 


0.999 


10 




Adopted model 


BPS CS 22956- 


-0094 


34 280 + 800 


5 63 + 2 


0.622 




2 ± 1 


SFIT 






34 100 ± 2000 


5.52 ± 0.2 








Ionization equilibrium 






34 000 


5.5 


0.699 


5 




Adopted model 


SB 21 




35 960 + 500 


5.4 ± 0.2 


0.997 




12 ± 2 


SFIT 






36 500 ± 1 500 


5.6 ± 0.2 








Ionization equilibrium 






36 000 


5.5 


0.999 


10 




Adopted model 






35 000 


5.4 








Huneeretal. (19811 


LB 1766 




36 340 ± 500 


5.19 ±0.1 


0.997 




20 ±3 


SFIT 






35 600 ± 2 100 


5.15 ±0.25 








Ionization equilibrium 






36 000 


5.0 


0.999 


10 




Adopted model 






38 000 ± 500 


5.5 ± 0.3 








Ahmad & Jefferv (2003) 






40 000 


6.3 








Lanz et al. (2004) 


BPS CS 29496- 


-0010 


39 150 ± 1000 


5.65 ±0.2 


0.996 




2±1 


SFIT 






39 770 ± 2300 


5.8 ± 0.2 








Ionization equilibrium 






40 000 


5.5 


0.999 


10 




Adopted model 


LB 3229 




40 000 ± 500 


5.15 ±0.2 


0.988 




8.5 ±2 


SFIT 






39 800 ± 1700 


5.34 ±0.4 








Ionization equilibrium 






40 000 


5.0 


0.989 


10 




Adopted model 



6) In all stars (four out of six) where it can be measured, neon ap- 
pears to be significantly overabundant. 

7) The Hel line fits are not uniformly satisfactory. In many cases, 
the observed line cores are substantially stronger than in the theo- 
retical profiles. Here we suspect possible non-LTE effects. In a few 
cases, particularly where the S/N ratio is low, one or both of the 
observed line wings lies below the theoretical profile. Here we sus- 
pect difficulties with the normalisation - which has always been as 
conservative as possible. 



4 EVOLUTIONARY STATUS OF HE-SDBS 

The evolutionary status of He-sdB and He-sdO stars has been dis- 
cussed for some three decades. The location of these stars on and 
above the helium main-sequence points to low-mass helium stars 
which are either currently in a helium-core burning phase (the he- 
lium main sequence), or are approaching or leaving such a phase. 

There has been a tendency to a nalyze and interp ret He-sdO 
l Stroer et a l. 2007; Napiwotzki 2008) and He-sdB stars dLanz et all 
|2004 Ahmad & Jeffery 120061) separately. This is partly due to the 
increasing importance of NLTE effects in sdO stars, and hence the 
adoption of different model atmospheres; a more self-consistent 
treatment of the two groups would be valuable. The apparent sepa- 
ration of the two groups in T c s — \ogg space (Fig. [3} may simply 
be a consequence of a significant range of T c g within a group of 
stars having comparable luminosities. Nevertheless, it is clear that 
extremely helium-rich subdwarfs mark out a clear locus in this dia- 
gram, and that they are more luminous than the extended horizontal 
branch where hydrogen-rich sdB stars are generally found. In addi- 
tion, there are a small number of cooler helium-rich stars, including 
V652 Her and BX Cir which may be closely linked to the helium- 
rich subdwarfs (He-sd's). 

The question of binarity is unresolved. The prototype He- 
sdB star PG 1544+4 8 8 is a binary consisting of two He-sdB stars 
jAhmad et aI1l2004l) . Only one other binary H e-sd, the double He- 
sdO star HE 0301-3039 jLisker et al.l |2004|) has been reported. 




Figure 3. Comparison of hot helium-rich subdwarfs and extreme helium 
stars with the evolution o f a helium white dwarf with Z = 0.001 following 
a late helium core flash | Miller Bertolami et al. 2008). The locations of nor- 
mal sdB's tedelmann et al.ll2003l) He-sdO's |Str6er et al.ll2007h . He-sdB 's 
(Ahmad & Je fferv 2003), EHe's IPandev et al.l2006l) and the stars discussed 
in this paper (N He-sdB) are identified separately (see key). For the latter, 
we adopt the SFIT results cited in Table 3. 



In such cases the extreme surface composition is a clear conse- 
quence of a close binary interaction, probably following a common- 
envelope phase, in which the entire hydrogen-rich envelope has 
been ejected. The question is then what differentiates the produc- 
tion of an He-r ich subdwarf from a c onventional H-rich sdB star in 
a close binary. Ijustham et al.l d2010j) propose a possible evolution 
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Table 3. Chemical abundances 



Star 


Ref 


H 


He 


C 


N 


O 


Ne 


LB 1766 




< 8.5 


11.54 


7.10 ±0.23 


8.29 ±0.25 


7.13 ± 0.30 


8.62 ±0.42 


SB 21 




< 8.5 


11.54 


6.73 ±0.18 


8.24 ±0.22 


7.25 ± 0.27 


8.49 ± 0.47 


BPS CS 29496-0010 




< 8.5 


11.54 


6.88 ±0.20 


8.48 ± 0.24 






BPS CS 22940-0009 




9.1 ±0.2 


11.54 


8.94 ±0.35 


8.46 ±0.22 


7.11 ± 0.34 


8.27 ±0.45 


LB 3229 




9.3 ±0.2 


11.54 


7.51 ±0.27 


8.66 ±0.28 




8.62 ±0.38 


BPS CS 22956-0094 




11.1 ±0.2 


11.40 


8.52 ±0.35 


8.33 ±0.24 






LB 1766 


1 


9.53 


11.53 


7.06 


8.77 






BXCir 


2 


8.1 


11.5 


9.02 


8.26 


8.04 




V652 Her 


3 


9.38 


11.54 


7.14 


8.93 


7.54 


8.38 


JL87 


4 


11.62 


11.26 


8.83 


8.77 


8.60 


8.31 


LS IV-14°116 


5 


11.95 


11.23 


8.47 


8.23 






Sun 


6.7 


12.00 


[10.93] 


8.52 


7.92 


8.83 


[8.08] 



Star 


Ref Mg 


Al 


Si 


S 


LB 1766 


7.17 ±0.09 


6.20 ±0.18 


7.03 ±0.18 


6.71 ±0.20 


SB 21 


7.24 ±0.10 


6.22 ± 0.10 


6.98 ± 0.26 


6.54 ±0.13 


BPS CS 29496-0010 


7.80 ±0.19 




7.05 ±0.23 


6.65 ±0.15 


BPS CS 22940-0009 


7.27 ±0.18 


6.12 ± 0.15 


7.23 ±0.24 


6.45 ±0.15 


LB 3229 


8.21 ±0.15 




7.38 ±0.22 




BPS CS 22956-0094 


7.22 ±0.13 




6.82 ±0.14 


6.96 ±0.10 


BXCir 


2 7.17 


6.04 


6.91 


6.67 


V652 Her 


3 7.76 


6.49 


7.49 


7.44 


JL87 


4 7.36 


6.28 


7.22 


6.88 


LS IV- 14° 116 


5 6.95 




6.63 




Sun 


6 7.58 


6.47 


7.55 


7.33 



References: 1. Lanz et al. (2004), 2. Drilling et al. 1 1998) 3. Jeffery et al. 1 1999) , 4. Ahmad et al. (2007), 5. Naslim et al." 
feOld) , 6. iGreves'se^r^auva' ll { 1998b . ^ Ibziembowskf l 1998) [the solar helium abundance is the asteroseismic value for the 
outer convection zone, the solar neon abundance is the meteoritic value; other solar abundances are for the solar photosphere] . 











: # He-sdE 

^He-sdB (this poper) 
1 »sdB 
: T sd0 
: A Hs-sdO 

,EHe 
I ^V652 Her 




/ ■ 

■ ■ 
■ 

■ y 




- jy 

: 




ZAMS / 





» **. • 

»i \^ HeMS 


* / :ahb 

• / 






1 








, 1 , , , 1 


f T T T T T T 1 




5.0 4.8 


4.6 4.4 

log Teff 


4.2 4.0 


3 



Figure 4. Comparison of hot helium-rich subdwarfs and extreme helium 
stars with the evolution of a Mi = 0.3 Mq helium white dwarf follow- 
ing the accretion of .2 Mq helium to represent a double He-WD merger 
I Saio & Jeffery 2000). Other symbols as for Fig.[3] 



following a double-core common-envelope phase involving inter- 
mediate mass stars. 

In the case of a close binary where envelope ejection exposes 



CNO-processed helium, the enhancement of nitrogen is easily ex- 
plained by the conversion, and hence depletion, of carbon and oxy- 
gen in the CNO-cycle. A similar abund ance pattern is exhi bited by 
the low-luminosity EHe star V652 Her jjeffervetal.lll999h . 

N-rich He-sd's are less easy to explain in the absence of a bi- 
nary companion. Equally problematic is C-enrichment in either the 
single- or binary-star cases, since it requires t he addition of carbon 
from 3a-burning. The low-i EHe star BX Cir jDrilling etalj 19981) 
provides a C-rich analogy to V652 Her. Neon is normally produced 
by a-captures onto 14 N. Since both are plentiful in CNO-processed 
helium, a high-temperature episode in the formation of the He-sd 
might naturally give rise to an overabundance of neon. 

Three evolutionary models have been proposed to address 
these questions for single He-sd's. 



4.1 The late hot flasher 

iBrown et alj feOOll) proposed that single star evolution with en- 
hanced mass loss close to the tip of the RGB will produce a star 
that suffers its helium-core flash late on the white dwarf cooling 
track. Since the flash occurs off-center and when the outer layers 
are compact, flash-driven convection leads to mixing of the rem- 
nant H-envelope with the helium core, and possibly also with some 
carbon from the He-flash itself. The star initially expands to be- 
come a yellow giant, and then contracts towards the He main se- 
quence as the helium-burning layers migrates to the center of the 
star. Subseque nt calculations have examine d a number of variants 
of this model jMiller Bertolami et al.|[200^> (Fig. [3). The result is 
either a N-rich or a C-rich He-sdB. 
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4.2 The double helium white dwarf merger 

The merger of two helium white dwarfs has been proposed var- 
iously to account for both normal and helium-rich hot subd- 
warfs dWebb"inklll984l ; llbenl 1 1990l : rTiitukov & Yungel'Sonlll99Cl : 
ISaio & Jeffervll2000l) . The progenitors are considered to be short- 
period systems from which most of the hydrogen and angular mo- 
mentum has been ejected during a common-envelope phase. The 
surviving mass of hydrogen is very small relative to the total mixed 
mass (i.e. tha t of the disrupted white d warf) mHe/mmked > 1.4 x 
10~ 4 /0.296 llben & Tutukovl d 19860 . where the hydrogen mass 
fract ion would probably in crease for lower mass white dwarfs. 

ISaio& Jeffervl (2000) investigated the evolution of a helium 
white dwarf which rapidly accretes helium, i.e. as a result of merger 
with another helium white dwarf. Following off-center helium igni- 
tion, the star expands to become a yellow giant, and then contracts 
as the helium shell burns inward through a series of mild flashes. 
The stable end-state of such an evolution would most probably be 
an He-sdB or He-sdO star. The surface layers of such a star should 
be dominated by the nitrogen-enriched hel ium from the disrupted 
helium white dwarf. ISaio & Jeffervl d2OO0l) found no evidence for 
surface carbon enrichment since each shell flash produces very lit- 
tle carbon and the subsequent flash-driven convection reaches the 
surface only after the first He-shell flash. The evolutionary track 
from one such model is shown in Fig. [4] ISaio & Jef ferv ( 2000) 
found that one such model would successfully account for the ob- 
served properties of the pulsating EHe star V652 Her, which must 
be in the shell-flashing phase. 

4.3 The helium white dwarf plus hot subdwarf merger 

ljustham et ai] {2010) have proposed a model in which a close bi- 
nary containing a post-sdB star and a helium white dwarf merge. 
The post-sdB star is essentially a 0.46 Mq hybrid white dwarf con- 
taining a small (« 0.3 Mq) carbon-oxygen core and a helium enve- 
lope. The addition of fresh helium reignites the helium shell, and re- 
turns the star close to the helium main-sequence with a helium-rich 
surface. Population synthesis calculations yie ld a T e a — log q di s- 
tribution similar to that of the He-sdO stars o f lStroer et al] (2007). 

4.4 Evolutionary status of He-sdB stars 

Comparison of the locu s of He-sdB stars analysed here and by 
Ahma d & Jeffervl d2003l) with the evolutionary calculations dis- 
cussed above (Figs.[3]and[4j suggests that only the double helium 
white dwarf merger model successfully accounts for the distribu- 
tion and surface composition of the N-rich He-sdBs. However, the 
evolution of both late hot flashers and white dwarf mergers is af- 
fected strongly by the metallicity, mass and envelope hydrogen con- 
tent. Further exploration of the parameter space and of population 
statistics will be necessary to explain the origins of these evolved 
stars. 

4.5 Surface-chemistry evolution in He-sdB stars 

In all of the above evolutionary models, the initial surface composi- 
tion might be assumed to be determined by the mean mass fractions 
obtained by combining several layers of stellar material. 

For example, the surface of the double helium white-dwarf 
merger would comprise any surviving hydrogen on the progenitor 
white dwarfs mixed with the helium core of the less massive com- 
ponent. The surface hydrogen to helium ratio (by mass) would then 
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Figure 5. Distribution of helium abundance versus T c ff and log g. The he- 
lium abundance logj/ = log(nn e /nn) is indicated by colour (or grey- 
scale) as shown in the key. The data shown include the He-sdB, He-sdO, 
sdB and sdO atars from Figs.[3]and|4] 



correspond roughly to the mass ratio of the hydrogen and helium 
layers in the latter, i.e. w 5 x 10~ 4 (see above, = 2 x 10~ 3 by 
number). Meanwhile the CNO abundances would lie somewhere 
between a fully CNO-cycled mixture (i.e. carbon and oxygen con- 
verted to nitrogen) and a primordial mixture (i.e. carbon, nitrogen 
and oxygen scaled to the iron abundance), depending on the he- 
lium to hydrogen fraction. Additional carbon might be present if 
3a burning occurs during the merger. 

The surface of a late-flasher will also be represented by 
CNO-processed helium, doped by whatever envelope hydrogen re- 
mained on the surface of the giant before hydrogen-burning was 
extinguished. This may be anywhere in the ranges rm/riHc ~ 
0.001 - 0.01 (shallow mixin g) or 10" 5 - 10" 6 (deep mixing) 
(Miller Be rtolami et al. 2008). Carbon-enrichment may occur if 
flash-driven convection can drive material to the surfce after 3a 
ignition, and if the helium-envelope mass is small compared with 
the amount of carbon available for mixing (Cassisi & Vinkll2003l: 
lLanz et al.ll2004l) . 

The question is then what happens to this mixture as the star 
contracts towards the helium main sequence. It is well known that 
in a radiative stellar atmosphere having a sufficiently high surface 
gravity, an imbalance between the relative radiative and gravita- 
tional forces on individual ions leads to chemical separation of 
atomic species by the process of diffusion. Diffusion in su bdwarf 
B sta rs, with log 5 > 5.5 and 25 000 < T cft /K < 40 000 faeberj 
1 19921) causes hydrogen to float and helium to sink so that observed 
helium abundances are typically in the range 0.001 < tihc/^h < 
0.01. 

At least two factors moderate the instantaneous conver- 
sion of a contracting He-sd into a normal H-rich sdB star. The 
first is that the diffusi on timescale is relatively long w 10 J y 
dUnglaub & Buesll200ll. no stellar wind). The evolution timescale 
essentially the thermal timescale for the enve lope and is m 10 6 y fo r 
a merger involving a O.2M0 white dwarf dSaio &Jeffervi r2000), 
or ss 10 4 y for a late flashe r with a 0.02Mq remnant envelope 
dMiller Bertolami et al.ll2008l) . 

The second factor is that a stellar win d acts to slow t he dif - 
fusion process to give a timescale ~ 10 y dUnglaubll2005l t2008). 
Since winds in hot subdwarfs are radiatively driven, they are lumi- 
nosity sensitive. Hence diffusion becomes more effective at low 
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luminosity (high gravity). Helium depletion in the photosphere 
will accelerate as a star contracts towards the zero-age horizontal 
branch. 

In the absence of winds, and if all He-sd's were formed in 
the same way with identical envelope masses, one might expect to 
see a helium abundance gradient along the observed sequence, or 
at least to see the helium abundance drop as a helium-rich sub- 
dwarf contracts across the "wind-line" - some critical luminos- 
ity below which diffusion becomes effective at reducing photo- 
spheric helium. The fact that the most hydrogen-rich "He-sds" {i.e. 
BPSCS 22956-0094, JL87, LSIV-14°116) are also those closest 
to the sdB domain might be important in this regard. 

To explore such an hypothesis, Fig. [5] shows the distribu- 
tion of hot subdwarfs as a function of helium abundance log y = 
log(riHc/?iH); the sample is exclusively that given in Figs.[3]and[4] 
These observations indicate a substantial majority of normal sdB 
stars with negligible helium, a significant number of extremely 
He-rich stars on the pre-subdwarf cooling track and a continuum 
of hot subdwarfs with relatively high gravities and intermediate 
helium abundances (—1 < logy < +1). In numerical terms, 
there are roughly 100 normal sd Bs for every ten helium-rich plus 
intermediate-helium subdwarfs ( Green et alj|l9861) . and eighteen 
helium-rich He-sdBs for nine intermediate He-sdBs (Fig. [5}; bet- 
ter statistics would be valuable. 

Hot subdwarfs with intermediate helium abundances might 
then represent stars in which diffusion has started to operate, but 
in which some helium remains visible. The ratio of intermediate 
helium subdwarfs to normal subdwarfs should then be given ap- 
proximately by the ratio of the diffusion timescale (~ 10 y) to 
the sdB nuclear timescale (~ 10 8 y) times the fraction of sdB stars 
formed through channels which involve a helium-rich progenitor 
(< 0.5), leading to a total of roughly one intermediate helium sub- 
dwarf in one hundred, slightly fewer than observed. 

One problem with this argument, and there are many, is that 
helium-rich sdB stars may simply become helium-rich sdO stars 
on the helium main sequence and may not evolve into helium-poor 
sdB stars. A more detailed examination of surface abundances of 
other species, including iron, in several helium-rich hot subdwarfs 
will be necessary. 

One observation, howe ver, is instructive. Among s t the helium- 
rich s ubdwarfs studied by lAhmad & Jeffervl d2006h . I Stroer et al. 
(2007) and ourselves, stars with intermediate helium abundances 
lie predominantly at the boundary between the He-poor and He-rich 
subdwarfs in the log g — T c g diagram. There are virtually no He- 
poor subdwarfs significantly above the horizontal branch (Fig.O, 
i.e. with T c ft > 30 000K and log g < 5.3. This observation is sup- 
ported by low-resol ution classification surveys dWinter et al 120061 ; 
iDrilling et al]|2003l) . Thus, if subdwarfs evolve onto either the he- 
lium main sequence or the extended horizontal branch by contract- 
ing from a more expanded configuration, then the only ones which 
are currently observed to be doing so are helium rich. 



5 CONCLUSION 

As part of an extended study of the surface abundances of ex- 
tremely helium-rich hot subdwarfs, high-resolution optical echelle 
spectra of the He-sdB stars LB 1766, SB 21, BPS CS 22940-0009, 
BPS CS 29496-001, BPS CS 22956-0094 and LB 3229 have been 
presented. Opacity-sampled line-blanketed model atmospheres 
have been used to derive atmospheric properties and surface abun- 
dances. 



All the stars analysed are moderately metal poor compared 
with the Sun ([Fe/H]«-0.5). LB 1766 and SB 21, BPS CS 29496- 
001 and LB 3229 are nitrogen-rich He-sdBs, while BPS CS 22940- 
0009 and BPS CS 22956-0094 are carbon-rich He-sdBs. The for- 
mer have a surface composition and L/M ratio comparable with 
the extreme helium star V652 Her, while the latter might be more 
directly compared with the extreme helium star BX Cir. 

The evolutionary status of He-sdB 's has been discussed in the 
context of i) close-binary star evolution, ii) a late helium flash in a 
post-RGB star, iii) the merger of two helium white dwarfs, and iv) 
the merger of a helium white dwarf with a post-sdB star. The sur- 
face composition and locus of single N-rich He-sdBs are currently 
best explained by the merger of two helium white dwarfs, although 
this may not necessarily be an unique solution. The merger of a he- 
lium white dwarf with a post-sdB white dwarf offers an interesting 
alternative. C-rich He-sdBs require further investigation; the origin 
of surface carbon is difficult to explain without mixing 3a products 
to the surface and so far, only the late flasher model seems capable 
of this. An over-abundance of neon requires further explanation. 

On the basis of any of these evolution tracks, the EHe stars 
V652 Her and BX Cir are likely to evolve to become He-sdB stars. 
He-sdB stars are likely to evolve to become He-sdO stars. 
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APPENDIX A: ONLINE MATERIAL 



Abundance analyses of helium-rich subluminous B stars 1 1 

Table Al. Oscillator strength (logg/), measured equivalent width (w\) and derived elemental abundance ti for each line measured in the six programme 
stars. 
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Notes: 1: empirical oscillator strength to match observed line: not used in mean. 2: blended with Nil: not used in mean, except for LB 3229 where Nil 
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Figure Al. The AAT spectrum of LB 1766 along with the best-fit model of T cff = 36000 K, logg = 5.0, n Hc = 0.999, vsini = 20 km s -1 , and 
vt = lOkms . Abundances are as in Table[3] 
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Figure A2. The AAT spectrum of SB21 along with the best-fit model of T off = 36 000 K, logg = 5.5, n Ho = 0.999, vsini = 12kms _1 , and 
vt = 10 km s~ 1 . Abundances are as in Table[3] 
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Figure A3. The AAT spectrum of LB 3229 along with the best-fit model of T cff = 40 000 K, log# = 5.0, n Hc = 0.989, vsini = 8.5kms _1 , and 
vt = lOkms -1 . Abundances are as in Table[3] 
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Figure A4. The AAT spectrum of BPS 22940-0009 along with the best-fit model of T e g = 34 000 K, log g = 4.5, n He = 0.989, V sin i = 4kms~ 1 , and 
vt = lOkms . Abundances are as in Table[3] 
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Figure A6. The AAT spectrum of BPS 29496-0010 along with the best-fit model of T e g = 40 000 K, log g = 5.5, n Hc = 0.999, V sin i = 2kms~ 1 , and 
vt = 10kms — x . Abundances are as in Table[3] 
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Table A2. Abundance errors 6ei due to errors in T c g and log g. 
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